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Abstract 

We present a new quantum complexity class, called MQ 2 , which is 
contained in AWPP. This class has a compact and simple mathemati- 
cal definition, involving only polynomial-time computable functions and 
a unitarity condition. It contains both Deutsch-Jozsa's and Shor's algo- 
rithm, while its relation to BQP is unknown. This shows that in the 
complexity class hierarchy, BQP is not an extraordinary isolated island, 
but has "siblings" which as well can solve prime-factorization. 

1 Introduction 

Quantum computing used to be a very popular discipline ranging from pure 
theoretical questions, concerning the complexity of the quantum polynomial- 
time class BQP, to practical concerns such as how to build a quantum computer. 
Nowadays, it seems a little bit that scientists are loosing their interest. Is it 
because all " easy" questions have been answered and what reamins is too hard or 
unintresting? We can find many unanswered questions considering complexity 
classes. For example, we know that BPPCBQPCAWPP. But questions whether 
is BQP equal to its "father" AWPP or its "son" BPP have not been answered up 
to now. Here, wc do not answer them cither. Instead, we introduce a nontrivial 
" brother" , which we call MQ 2 . Surprisingly, this brother can also factorize long 
integers in polynomial time. Moreover, it has a very compact mathematical 
definition, which does not explicitly involve any physics. 

The paper is organized as follows. In section 2, we introduce the neccessary 
notation and definitions. In section 3, we briefly review classical polynomial- 
time classes definition. In the following section, we define the class MQ 2 it- 
self. Fifth sections demontrates two quantum algorithms, Shor's and Deutsch- 
Jozsa's, in class MQ 2 . 



2 Definitions and notation 

In this text, we frequently encounter matrices. All matrices here will have square 
shape. We will denote the element of a matrix M in i— th row and j—th column 
as (i|M|j), in accordance with the usual notation used in quantum computing, 
because we think it makes the text more readable than writing Mij. The range 
of the indices will be 0..n — 1 where n is the number of columns or rows. 

We will define our new complexity class with use of matrices. To make the 
class uniform, we will want that all the matrices are constructed by the same 
algorithm. We formalize this in the following definition. 
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Definition 1 (Poly-computable matrix family). A sequence of matrices 
T11T2, ■ ■ ■ is called a poly- computable matrix family if there exists a function 
/:NxNxi-»l computable in time polynomial in the length of all the 
arguments such that 

Vi,i,n : f(i,j,n) = (j\T n \i). 

Sometimes, we will drop the index n when it will be clear from the context. 

To show how the classical polynomial-time classes definitions correspond to 
our definition, we will start with a Turing machine and express its transition 
function as a transition matrix. A transition matrix is in fact a linear operator 
defined in a vector space spanned by configurations playing role of base vectors. 

Definition 2 (Configuration of a Turing machine). A configuration of a 
Turing machine is an ordered triple consisting of: 

-the contents of the tape 

-the current state 

-the position of the head 1 

We emphasize here that the configuration as defined above contains also 
the content of the tape, which is not true for configurations as defined else- 
where. Without loss of generality, we will further assume that configurations 
are indexed and denoted by their indices. A special position among the con- 
figurations has the family of initial configuration I(x), which we allow to be 
dependent on x, but require to be computable in polynomial time. Again, we 
will drop the index n when it will be clear which member of the family we mean. 
We will also need to be able to recognize accepting configurations. For this pur- 
pose, we will have a function a(x,c), computable in polynomial time, where 
the first argument will be the input for the algorithm and c is a configuration. 
The function a(x, c) will return 1 iff the configuration c is accepting (possibly 
depending on x) and otherwise. 

Now we are ready to jump to the notion of transition matrix. If a configu- 
ration ci leads to another configuration C2 in the next step with probability p, 
there is p on the position (c2|T|ci), otherwise there is zero. Because the tapes 
are of unbounded size, so is the matrix. However, if we know that the time 
complexity of a Turing machine is T{n), we may for fixed n have a finite matrix 
cutting the tapes at the distance T(n) from the initial position on both sides. 
The size of the matrix for inputs of length n is then 2°( T ("» x 2°( T (™». For 
a probabilistic Turing machine, the transition matrix is stochastic, e.g. every 
row sums up to 1 . Transition matrices naturally form a poly-computable matrix 
family, since for each pair c±, C2, the probability of going from one to another 
can be read from the description of the underlying probabilistic Turing machine, 
which is a finite object 2 . 

3 Traditional complexity classes 

We will now briefly review classical complexity classes definitions. The common 
definitions of P, BPP, NP and PP involve a probabilistic Turing machine and 
look at the accepting probability for each input 3 . 

1 We assume without loss of generality that the machine has only one tape. 

2 It should be pointed out that for most poly-computable stochastic matrix families, no 
corresponding probabilistic Turing machine exists. The reason is that each probabilistic Turing 
machine has finite description of a bounded size, independent on the length of input, while in 
our definition we allowed function / to arbitrarily depend on n. 

3 Obviously, the class P can be viewed as a special case of a probabilistic class, with prob- 
abilities either one or zero. 
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One step of a probabilistic Turing machine corresponds to multiplying the 
transition matrix with a vector representing the current configuration. Thus, 
instead of saying "the probability of accepting on a configuration I(x) after S 
steps is p", we may eqvivalently say "E C :a(c a: )=i( c l^ ,S |^( x )) = P"- We will use 
this observation in the following definitions. 

Definition 3 (Polynomial time classes in matrix notation). A language L 
is in class C if there exists a polynomial pin) and a probabilistic Turing machine 
M with transition matrix family Ti and functions I(x), a{x,c) computable in 
polynomial time, such that for all n and for all x of length n: 
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where exactly one of the columns applies. 



In the same manner, we may define the quantum class BQP: 

Definition 4 (BQP). A language L is in class BQP if there exists a polynomial 
p(n) and a quantum Turing machine M with transition matrix family Ti of 
unitary matrices and functions I{x), a(x,c) computable in polynomial time, 
such that 



For x G L : 
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We emphasize here that there are exactly two points in which Definition 0] 
and the Definition of BPP in Definition [3] differ: First, in Definition |21 we have 
stochastic matrices while in Definition 0] we have unitary matrices. Second, in 
Definition |21 we are looking at the value of ~^ c . a ( c x ) = i{c\T plyn ^\I{x)) , while in 
Definition 01 we look at the square norm of this value. However, the latter can 
be avoided, since we can eqvivalently use the square norm in the Definition |21 of 
BPP. Thus, the only remaining difference between the two classes is the type of 
matrices used. 

Now we will define the class MQ 2 itself. 



4 The definition 

We alter the Definition 0] of class BQP to get a new class MQ 2 . At first, we 
will not use transition matrices, but instead a unitary, poly-computable matrix 
family. That is a weaker requirement. In turn, we will be more strict in the 
number of matrices allowed - we will only use two copies of the matrix. 
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Definition 5 (MQ 2 ). A language L is in class MQ 2 iff there exists a uni- 
tary, poly- computable matrix family T, A poly- computable vector family I{x) 
and function a(cc, c) computable in polynomial time such that 4 



ForxGL, E c: „( c ,*)=l (c\T 2 \I(x)) 



> 



ForxiL, £c:a(c,*)=i (c\T 2 \I(x)) 



< 



The resulting hierarchy is visualised in Fig ureffl The inclusion MQ 2 CAWPP 
was shown in 0]. 
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Figure 1: Hierarchy of classes including MQ 2 . For each pair connected by a line, the 
class that stays upper contains the lower one. 



5 Expressing quantum algorithms 

In this section, we will demonstrate how two famous quantum algorithms fit 
into the class MQ 2 . Here we will only present the main ideas of the proofs. Full 
proofs can be found in 

At first, we show that class MQ 2 captures Deutsch-Jozsa's problem. For 
this problem, see In the Deutsch-Jozsa problem, a quantum oracle is used. 
That is a diagonal quantum gate, or in other words a diagonal unitary matrix, 
realizing the transformation x —* {— l)^ x \ Here, our matrix will be a product 
of a poly-computable matrix and this oracle. The result if then poly-computable 
too. 

Theorem 6. The class MQ 2 solves the Deutsch-Jozsa problem. 

4 The numbers 2 and i in the definition can be amplified in the same way as in the case 
of BPP and BQP. Proven in g| 
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Figure 2: The trick used to Et Deutsch-Jozsa's algorithm into the MQ 2 class. In a), 
there is the original setup, in b) there is the one we use, which gives the same result 
up to a global phase. 



Proof sketch. We will mimic the circuit from Deutsch-Jozsa's algorithm (see 
Figure El a)) by two copies of a poly-computable matrix family T (see Figure |2 
b)). For a fixed n, the matrix T will be a product of the oracle and H n . We 
may add another matrix for the / function to the front, since it will only add 
the number (— to the global phase and thus will not change the result. 
Formally, we define a matrix T as 



(y\T\x) °= {-l) f ^(y\H n \x) = _l = (_l)/(»)(_l)E 4 »*wniod2 

v2™ 



which is obviously computable in polytime and unitary. We define Cj(x) = 0" 
and ca(x) = 0™ for x of length n. Then we have 



|(0"|T 2 |0") = |^(0"|T|fc)(fc|T|0" 

k 

' . /on v ' 



1 2 

111. .,12 1 _^ ^ ^/(0") ^ i "i fe i m ° d 2 



f(k) 



If the function is constant, then the sum l)^^-* equals ±2™ and the prob- 

ability |(0™|T 2 |0 n )| 2 equals 1. If the function is balanced, both the sum and the 
probability is 0. 

□ 

Now we will show how to implement the famous Shor's algorithm |S] in class 
MQ 2 : 

Theorem 7. The class MQ 2 solves the factoring problem. More precisely, 
there exists a constant k such that given numbers x and N>k as in the Shor's 
algorithm, the language 

L = {(N, i)\x a mod N has a period r whose i-th bit is 1} 



is in MQ 2 . 
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Figure 3: The trick used to tit the Shor's algorithm into MQ 2 class. In a), there is 
the original setup, in b) there is the one we use, which gives the same result up to a 
global phase. 



Proof sketch. We will use the same notation as in the original paper by Shor [5] . 
N is the number to factorize and N — P1P2 where both p\ pi are primes and are 
different from each other. We then arbitrarily choose an x coprime to N. The 
pair (x, N) is the input of the algorithm. The goal is to find the smallest r ^ 
such that x r mod N = 1. This r is called a period of x. We choose a number q 
such that q is a power of 2 and q > 2 2 ^ l ° 92N \ This number will, together with 
the length of x and N, determine the size of the matrix. 

In the original setup, see Figure|3]a), we have three different transformations, 
namely: Hadamard transformation, a transformation computing power modulo 
N, and the DFT. One can notice that the effect on 0" of the DFT and of the 
Hadamard transform is the same, so we suffice with only two different transfor- 
mations. Furthermore, the matrices realizing them are poly-computable: 

(x',N',a',i'\DFT\x,N,a,i) °= ^-S x ^S N , N ,S^e 2 f Laa ' 

(a/, N', a', i'\MOD\x, N, a, i) D = F 5 XtX i5N,N>8 a ,a'8v+i,xa mod n 



One may also simply check that both the matrices are unitary. Their product 
T D = F DFT ■ M OD is thus unitary too. Its elements read 

(x ,N ,a ,i \T\x,N,a,i) = — — 5 x . x ,8 N . N >e * a a 5 i+i > tX a mod N 

and are again clearly functions computable in polytime. Applying two times the 
matrix T, we apply an extra MOD transformation comparing to the original 
setup. Nevertheless, on 0", this transform acts as identity. Any classical post- 
processing, as in the original algorithm, can be incorporated into the function 
a{x, c). 

□ 



6 Conclusion 

We saw a complexity class which has a compact purely mathematical definition. 
In order to describe an algorithm, we suffice with three polynomial-time com- 
putable functions taking bit-strings as arguments: f(i,j,n), I(x), and a(x,c). 
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Even quite a complex algorithm, as the Shor's certainly is, can be described on 
three lines. Further, the class MQ 2 shows that BQP is not the only possible 
class, lying in between BPP and AWPP, and not being trivially equal to either 
of them, which can do factorization and exponential speedup with oracles as in 
Deutsch-Jozsa's algorithm. 
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